N92-22 085 

TOWARDS A WAVE THEORY OF CHARGED BEAM TRANSPORT: 

A COLLECTION OF THOUGHTS 


G. Dattoli, C. Mari, and A. Torre 
ENEA - Area INN, Dip. Sviluppo Tecnologie di Punta, 
p.o. Box 65 - 00044 Frascati, Rome, Italy 


ABSTRACT 


In this paper we formulate in a rigorous way a wave theory of charged 
beam linear transport. The Wigner distribution function is introduced 

and provides the link with classical mechanics Finally *** ?!!" Tinllr 

equation is shown to coincide with the Liouville equation for the 1 

transport. 


INTRODUCTION 


A formal ''quantum" theory of charged beam transport has been recently 
proposed. (Ref. 1) Within such a context the possibility of viewing the 

beam emittance as a kind of quantization constant has J?®®" ^witt^the 
on the basis of some strong conceptual similarities existing with the 

so-called "quantum" theory of light rays. (Ref. 2) 

The proposed quantization procedure (described in some detail below) 
cannot be thought as fully satisfactory and in particular the role played 
by the "beam wave function" (b.w.f.) and its link with classical dynamics 
have not ben thoroughly investigated and clarified. 

In this paper we develop, in a more rigorous way, the fundamental 
steps towards a "quantum" theory of beam transport through linear elements^ 
showing the relation with classical mechanics and discussing the con 
nections with the Liouville equation which describes the time evolution 
of non-interacting classical ensembles. We will show that classical 
dynamics and the Courant-Snyder theory (Ref. 3) can be recovered by using 
the formalism of the Wigner phase-space function, (Ref . 4 _therefore 
can be seen as the natural framework to study the evolution of charged 
beam transport through linear magnetic systems. 

The paper is organized as follows. In Sect. 1 we review some aspects 
of symplectic mechanics (Ref. 5) and introduce some definitions which 
will be useful in Sect. 2 where the "quantum" theory of charged linear 
transport will be completely developed in terms of the Wigner phase space 
function, 'swne comments and final remarks on the -quantum" theory of 
nonlinear charged transport conclude this paper. 
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1. 


SYMPLECTIC MECHANICS OP QUADRATIC HAMILTONIANS 


The most general time dependent quadratic Hamiltonian in one degree 
of freedom can be written as 


w = ^ a (0q 2 + ^(0p 2 


c(t)qp 


or in a matrix form as 

1 T . 

H = -z T Hz 

2 - 


where 



H, 


a(0 c(t)\ 

c(0 6(07 


and the superscript " T " denotes transpose. 


( 1 . 1 ) 


( 1 -2) 


(1.3) 


The equations of motion for the vector ^ are obtained from the Poisson 
Brackets (P.B.) with the Hamiltonian (1.1) . The relevant rules are easily 
derived. It is quite straightforward to realize that 


{z,z 7 } = S = 


0 

- 1 



( 1 -4a) 


and that 

{z,z t Bz} = 2SBz * (1.4b) 


where S' is the unit symplectic matrix in two dimensions. According to 
the above rules the equations of motion of ^ are easily written, namely, 

z={z,z t Hz) = SHz (1.5) 

and therefore immediately integrated, thus yielding 

z(O-0(Oz(O) (1.6) 

where 


Note that (2 r ,z} = 0and therefore the above product cannot be, strictly 
speaking, considered a conventional P.B. In the case that the second 
term in the {,} -product is a scalar function [as in Eq. (1.4b)] then 
the {,} is a conventional P.B. 
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(1.7) 


0(t) = { exp 


S fl (t')dt' 
J o 


and (,K denotes time ordering for the classical evolution operator which 
is necessary when the commutator [SH(t),SH(L')] is different from zero. 

The above formalism is particularly useful to introduce the 
fluctuation tensor and its dynamical features, namely, let us pose 


£(Q =< z(t)z T (t) > -< z>< z T > 

where <,> denotes ensemble average. 

In a matrix form we recover the usual expression 


( 1 . 8 ) 


>.(0 = 


o. 


PQ 


0 

a 


2 ' 
PQ 
2 

P . 


(1.9) 


The equation of motion of Z follows from Eq. (1*5) and reads 


±yxt) = vi: + lv 

dt 


( 1 . 10 ) 


where V = S If , and can be immediately integrated in terms of the above 
evolution operator 0 thus getting 


£(f) = 0 (t) 


( 1 . 11 ) 


Within this formalism it is easy to recover the well-known quadratic 
invariant of Courant-Snyder , namely, let us write the most general 
quadratic (time dependent) expression in ^ 


I=z ! T(t)z 


( 1 . 12 ) 


where 


7(0 = 


Y a 
a 


— = —+ {/ ,// } = 0 
dt dt 


(1.13) 


Then we impose that / is a "time-dependent” invariant, i.e., 


(1.14) 


thus obtaining the following equations specifying a,p,y 
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(1.15) 


/p = -26a + 2c(3 P(0) = p o 

/ Y = 2aa - 2cy y(0) = y 0 
\^a = ap-pY a(0) = a 0 

from which it is easy to verify that 
PY = 1 + Qt 2 (1.16) 

The Courant-Snyder invariant reads explicitly as 

/ = Y<7 2 + 2apg + pp 2 (1.17) 

and it is widely used in the theory of linear transport. 

2. TOWARDS A QUANTIZED THEORY 


In the previous section we introduced the necessary background to 
derive a kind of uncertainty principle in the theory of linear transport. 
From Eq. (1.11) we get the relevant result 

detl(t) = det£(0) ' (2.1) 

i.e., det£ is an invariant quantity and the emittance A of the system 

^(0-[<<7 2 ><p 2 >-<<7P> 2 ]‘ /2 (2.2) 

is preserved in time. 

Furthermore this means that 

o q o p >A 


which represents a kind of uncertainty principle in the canonical variables 
q,p and can be used as the starting point of our quantization procedure. 

rules are simple, the beam momentum is replaced by an operator 
specified by 

= ~Ard_ 

P ~~iYq (2.3) 


where -^Pis the beam reduced emittance and the following rule of commutation 
is also assumed: 


[p,g] = -iyp 


(2.4) 
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A "Hamiltonian" operator is finally associated to the longitudinal 
coordinate of propagation s 


H = ivh- 


ds 


(2.5) 


so that the "Schrodinger" equation for a beam passing through a quadrupole 
of strength /c(s) reads 


d 


ds 


vi 2 a 2 

2 dq 2 


^k(s)q 2 


V(t7.s) 


( 2 . 6 ) 


where \ p(g,s) denotes the beam wave function (b.w.f.). Clearly the b.w.f. 
ip(g,s) must be related in some way to the "classical" beam distribution 
p (g.p.s) satisfying the Liouville equation for the time evolution of an 
ensemble of single-particle systems. The link is not obvious since op 
depends only on q and eventually on s while p is a function of q , p and 
s.The answer is given by the Wigner distribution function which is defined 
as follows: 


W(q,p,s)= 



g + -y,s hp < 7 --y,s e 


ipy / 


(2.7) 


and satisfies the Von-Neumann equation for a generic potential 1 7(g) (Ref. 

6 ) 


— + p — 

ds H dq 


W(q, p , s) =-yK 


v\ q ^±- p \-v 


2 dp 


l/(<7,p,s) 


( 2 . 8 ) 


In the above-considered case of propagation through a quadrupole of 
strength k{s ) i. e. , for an elastic potential 17(g) = ^(s)g 2 , the Von-Nuemann 


equation reduces to 


— l/(g,p;s) = - 
ds 


d d 

p — -k(s)q — 

L dq d p 


W(q,p\s) 


(2.9) 


which is equivalent to the Liouville classical equation. In this simple 
case we can identify the classical distribution density p with the Wigner 
distribution function W thus providing a completely consistent quan- 
tization scheme. Within this framework the physically meaningful quantity 
is not | ip(g , s) | 2 dg but W (g , p \s)dqd p which ensures the consistency of 
our procedure. 

As a simple exercise it is straightforward to see that 
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W(q,p;s) = 


I 

(2n)(<rt/2) eXP 


\ v(s)q 2 + 2a(s)qp->-ft(s)p 2 
2 (A/2) 


( 2 . 10 ) 


is a solution of (2.9) if and only if the "time-dependent" parameters 
(a,|3,y) satisfy the following system of differential equations: 

(p' = -2a 

<a' = *(s)p-Y (2.11) 

V y' = 2/c(s)a 


which are the well-known equations defining the evolution of the Twiss 
parameters in quadrupole lenses [see Eq. (1.15)]. 


3 . CONCLUDING REMARKS 

The extension of the developed theory to the case of nonlinear 
transport of charged beams is not straightforward and the problems involved 
can be illustrated in a simple example. 

A sextupole term introduces, in the single-particle Hamiltonian, a 
contribution of the type 

^(Q.s) = ^4pg 3 (3.1) 


Inserting the above potential in the Von-Neumann equation (2.8) we easily 
get the following evolution equation for the Wigner distribution function: 


-W(q,p;s) = - 


i r - i 

L *p J 


W(q,p\s) 


L 




12 


dp 


-W(q,p-,s) 


(3.2) 


where L denotes the Liouville operator associated to the potential (3.1) 
and given by L = - p^ + X( 7 2 ^. The extra term in Eq. (3.2) is a purely 

"quantum" contribution and it is not present in the classical Liouville 
equation for p (q , p\s) . From this point of view W and p cannot be identified 
and since they coincide in the limit vt -* 0, p may be viewed as the "classical" 
counterpart of the Wigner distribution function W . 
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